
%不同情况，需要更改调用的函数
[t,w]=ode23(@ifun,[0:0.1:20],[0,0.01,0,6,0,0.01],[ ]);
plot3(w(:,1),w(:,3),w(:,5));%绘制轨迹
grid on
title('b\neq0,e\neq0')   %标题每次改一下